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Abstract 


Starting from the finite volume form factors of local operators, we show how and un¬ 
der which hypothesis the c - 1 free boson conformal held theory in two-dimensions 
emerges as an effective theory governing the large-distance regime of multi-point cor¬ 
relation functions in a large class of one dimensional massless quantum Hamiltoni¬ 
ans. In our approach, in the large-distance critical regime, the local operators of the 
initial model are represented by well suited vertex operators associated to the free bo¬ 
son model. This provides an effective held theoretic description of the large distance 
behaviour of correlation functions in ID quantum critical models. We develop this de¬ 
scription starting from the hrst principles and directly at the microscopic level, namely 
in terms of the properties of the hnite volume matrix elements of local operators. 


Introduction 

It widely believed that the spectrum of a quantum Hamiltonian is intimately related to certain overall properties 
of the large-distance asymptotic behaviour of its correlation functions. If, in the infinite volume L limit, a given 
model’s spectrum exhibits a gap between sl finitely degenerated ground state and the tower of excited states above 
it, than the correlation functions are expected to decay exponentially fast in the distance of separation between 
the various operators involved in the correlator. This behaviour changes for gapless models, viz- those whose 
ground state, in the L -t-cx? limit, becomes directly connected to the continuum of excited states above it. For 
those models, the correlation functions are expected to decay algebraically in the distance of separation. The 
powers of the distance which arise in this algebraic decay are called critical exponents. Models having the same 
values of their critical exponents are said to belong to the same universality class. In fact, it is believed that the 
features determining a given universality class are rather sparse in the sense that solely certain overall symmetries 
of the model should fix if ifTTl . Unforfunafely, fhe lack of any exacfly solvable and fruly inferacfing many-body 
quanfum crifical model in dimensions higher fhan one did nol allow, so far, for a direcf and explicif check of fhese 
properfies. However, extensive numerical dafa and experimenfal resulfs do speak in favour of such an interplay 
1221 . 
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The situation improves drastically in the case of one spatial dimensiorEl where a plethora of exactly solvable 
models arises: the Luttinger model (see GOl and references therein), 1 + 1 dimensional conformal field theories 
IS [13 or models that are solvable by one of the versions of the Bethe Ansatz |ih| such as the celebrated XXZ spin- 
1/2 chain or the ID Bose gas at arbitrary repulsive coupling. The abundance of exact solutions turns these one¬ 
dimensional models into a laboratory allowing one to test the universality principle in many concrete situations. 

It is worth reminding that, in all these models, two scales coexist: a microscopic scale 6 which is related to the 
lattice spacing or, more generally proportional to the inverse of the Fermi momentum, and a macroscopic scale L 
corresponding to the volume or number of sites of the model. Scale invariance can only be realized at distances Ax 
between local operators that range between these two scales and but far from them. This imposes that the ratio of 
the two scales L/d is very large in such a way that there exist a whole range of distances verifying 6 « Ax « L. 
For Ax close to either the microscopic or the macroscopic scales the scaling properties of the correlation functions 
get modified in drastic ways. 

The very sfrucfure of a conformal field fheory imposes a simple form for fhe fransformafion of ifs operators 
under scaling. In ifs furn, fhis imposes fhaf fhe correlation funclion exhibif an algebraic in fhe disfance pre- 
factor. In particular, fwo-poinf functions are purely-algebraic. The exponenfs driving fhe power-law behaviour of 
a correlator are fhen construcfed in ferms of fhe scaling dimensions of fhe operators fhaf are being averaged. The 
dafa issuing from conformal field fheories can fhus be used so as to provide one wifh quife explicif predictions for 
fhe large-disfance decay of correlators in massless quanfum models in one spafial dimension. In facf, fhere, one 
can be slighfly more precise in respect to the two pillars on which these predictions build 

• in the critical regime 6 « Ax « L, correlation functions should exhibit conformal invariance as argued by 
Polyakov 1431 . This suggests that the leading contribution to the long-distance asymptotics should be repro¬ 
duced by correlation functions of appropriate operators in a two-dimensional conformal field fheory. Sfill, 
the identification of which conformal field fheory is fo be used and which are fhe "appropriafe" operafors is 
leff open af fhis sfage. 

• The 1/L corrections fo fhe model’s ground and low-lying excifed sfafe’s energies confain fhe informafion 
on fhe cenfral charge of fhe conformal field fheory describing fhe asymptotics and fhe scaling dimensions 
of fhe fower of operafors describing fhe correlafion funcfions as argued by Cardy 171191. 

• The choice of "appropriafe" operafors, on fhe conformal field fheory side, is done by advocafing fhaf fhese 
should inherif fhe symmefry fhaf is safisfied by fhe correlafion function one sfarfs wifh. 

To rephrase fhe above, Polyakov’s argumenf jusfifies why a conformal field fheory should emerge as an effec¬ 
tive large-disfance fheory in fhe domain 6 « Ax « L of fhe model while Cardy’s observafion permifs one fo 
exfracf, from fhe knowledge of fhe sfrucfure of a model’s excifafions, fhe quanfifies which would characferise fhe 
effecfive conformal field fheory describing fhe model’s large-disfance regime. 

A similar line of fhough is followed by exploifing fhe Luffinger liquid model as a fool for providing fhe critical 
exponenfs. Is is argued fhaf a model belongs fo fhe Luffinger liquid universalify class if if has fhe same form of fhe 
low-lying excifafions above ifs ground sfafe. The paramefers describing fhese excifafions fix fhe Luffinger liquid 
model fhaf is pertinenf for fhe model of inferesf. The critical exponenfs of fhe model one sfarfs wifh are fhen 
deduced from fhe ones compufed explicifly for fhe Luffinger liquid model ifTSlfT^IT^ . 

Independenfly of which of fhe above mefhods one chooses fo employ, one needs fo access fo fhe 1 /L correc¬ 
tions fo fhe energy specfrum of a model’s Hamiltonian so as fo be able fo carry ouf effecfive predicfions. The 
extracfion of such corrections from fhe specfra of various quanfum infegrable models was inifiafed in fhe series of 
works ifTOl fin [29] [30l |3T1 and led fo fhe idenfificafion of fhe cenfral charge and scaling dimensions, for numerous 

tin the language of classical statistical mechanics, this corresponds to the case of a two-dimensional model 1481 . 
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quantum integrable models. It is worthy mentioning the work lITSl where the 1 /L corrections to the free energy 
of certain non-integrable perturbations of the two-dimensional Ising model were obtained, on rigorous grounds 
through constructive field theory methods. By using Cardy’s form of the 1/L corrections, this work demonstrated 
that, for sufficiently small perturbations, the model still has the same central charge 1 /2 as the Ising model. 

Although effective in the sense that producing explicit answers, the above techniques are more of a list of 
prescriptions than a well argued from the first principles line of though that joins, argument after argument, 
the structure of a given microscopic model with the resurgence of an effective field theoretic description in the 
large-distance regime. Several attempts have been made in the literature to bring some elements of rigour or, at 
least, some ab inicio justification of the principle. In a series of works Il42ll44]| . by using the tools developed in 
Boll . Polyakov argued the arisal of a universal behaviour on the basis of perturbative field theoretic calculations. 
Although constituting an important step forward, his reasonings did not allow for any control on the magnitude 
of the contributions that he dropped from his calculations on the basis of some hand-waving arguments, without 
mentioning the problems inherent to lacks of convergence in perturbative handling in quantum field theory. A 
rigorous approach allowing one to prove the power-law decay of certain two-point functions in interacting one¬ 
dimensional fermion models that are a sufficiently small perturbation of a free fermion model was hrst proposed 
by Pinson and Spencer HTTl . further developed by Mastropietro OTl and then generalised in OBI to the case of 
multi-point energy-energy correlation functions in certain non-integrable perturbations of the two-dimensional 
Ising model. An adaptation of this approach |5l also allowed to establish the universality of certain properties 
of Luttinger liquid type for a class of sufficently small perturbations of a free fermion model. The method relies 
on the possibility to provide a rigorous construction of the path integral for such models in hnite volume and then 
study its inhnite volume and scaling limits through rigorous renormalisation group methods. Finally, we should 
also mention the approach for proving the universality of certain types of bi-dimensional Ising lattcies -which are 
related to one-dimensional massless quantum Hamiltonians by means of the correspondence proposed in ll48l - 
that was developed by Smirnov in t47l . 

The present paper introduces a convenient microscopic description for the spectrum, space of states and matrix 
elements of local operators in a class of one-dimensional massless quantum models. Our setting allows us to 
construct, in the large-distance regime 6 « Ax << L, a one-to-one map (up to higher order corrections in 
the distance) between the relevant to the large-distance regime sub-space of the model’s Hilbert space and the 
free boson Hilbert space. The local operators of the original model are then represented in terms of a collection 
of vertex operators acting in the free boson Hilbert space. The setting we introduce, as shown in the works 
122 [H [3a, is clearly verihed for a large-class of quantum integrable models, the XXZ spin-1/2 chain being the 
most prominent example. However, we do trust that the overall hypothesis that we lay down for the structure of 
the model’s observable is quite universal and at least encompasses several instances of one dimensional quantum 
liquids ifT^ . This is supported, e.g. by perturbative calculations around a free model 111 121 El- Furthermore all 
objects that we handle are standard within the phenomenological approach to interacting Fermi systems lfT^l38ll . 

Going slightly deeper into the details of the physics’ jargon, the main assumptions on which our setting builds 
are that 

• in the large-volume limit, the relevant part to the critical regime of the model’s spectrum is purely con¬ 
structed in terms of particle-hole excitation^; 

• the form factors -expectation values of local operators- taken between two states realised in terms of exci¬ 
tations that are ah located in an immediate vicinity of the Fermi zone- admit a structure descending from a 

^In the following, for simplicity, we shall assume a purely particle-hole excitation spectrum. This allows one to lighten the discussions. 
We stress that this assumption is not a limitation in that, for instance, the bound states are expected to solely produce exponentially small 
corrections in the distance. 
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local repulsion principle between the interacting momenta of the particle and hole building up the excited 
state. 

The main result of the paper can be phrased as follow. Let 0\{x\),... ,Or{Xr) be local operators located at at 
positions x\,...,Xr and associated with some one-dimensional quantum model in finite (but large) volume L and 
having the Hilbert space bphys- The operators Os is assumed to induce solely transitions between spin sectors 
differing by some fixed infeger Os- Then, in fhe large 5 \xa - Xb\ L limif, fhe r-poinf ground-lo-ground sfafe 
expecfafion value of fhese operafors satisfies 

(0.1) 

' ^ Ophys ' ^ Oetf 

where fhe - sign is fo be understood as an equalify up fo fhe firsf leading correction arising in each oscillating 
harmonics in fhe spacing difference. The operafor ffs, ^ - ■ ,r, appearing in fhe rhs of (10.11) acf in fhe Hilbert 

space heff associated to fhe free boson model while fhe expecfafion value on fhe Ihs of (10.11) is faken in fhe Hilbert 
space bphys of fhe initial model; fhey are builf up from fhe free boson vertex operafors Yiv, k', co): 

^sicos) = 2]r.(a)- 

KeZ 

In fhis formula, 

• is the variable attached to the left (-) or right (-I-) end of the Fermi zone which parametrises 
the position of the operator on the free boson model side. The factors cr* are certain re-scaling coefficients 
which take into account the variation of the dressed momentum at the endpoints of the Fermi zone. The are 
expressed in terms of the physical observables of the model, cf. (12.221) : 

• the indices L and R refer to the left or right copy of a free boson Hilbert space f) on which the vertex operator 
acts, viz. heft = hi ® hi?; 

• y {v, K \ io) is the vertex operator, cf. (11.271) . which connects states belonging to sectors differing by fhe 
charge a-; 

• Vi(±< 7 ) are fhe values of fhe so-called shiff funcfion relatively fo fhe ground sfafe, cf. (12.291) . and evaluated 
af fhe lefl (-q) or righf (q) end of fhe Fermi zone. We do sfress fhaf if is a funcfion of fhe charge ic and of fhe 
operator’s spin o^. 

• pivsiq) - Os + k) and p{Vs{-q) + k) is fhe fower of scaling dimensions associated wifh fhe operafor Os', 

• yK{Os), A" € Z, is fhe fower of fypical form facfors of fhe operafor Os faken befween fhe ground sfafe and fhe 
lowesf-lying excited sfafe belonging fo fhe A:-Umklapp secfor. 

We refer fhe reader fo Sub-Section 13. 3. II where explicif examples for fhe above-listed quanfifies are given in fhe 
case of fhe fundamenlal local operafors for fhe XXZ chain. In Sub-Secfion l3.3.2l we show how fhe strucfure of fhe 
Luffinger-liquid critical exponenfs can be recovered wifhin our formalism. 

Whaf equation (10.11) slates is lhal, in fhe weak sense -i.e. for expecfafion values- and for large disfances of 
separafion all fhe essence of fhe physical operator Oa{Xa) acfing in fhe Hilbert space bphys is grasped by fhe effecfive 
operafor ffsi^^s) acfing in fhe free boson model’s Hilbert space beff- Formula (10.21) provides an effecfive map, in 
fhe crilical regime, befween local operafors of fhe physical models and vertex operafors in fhe corresponding 
conformal field fheory. This map is completely defermined by fhe fypical form facfors fifOs), a- € Z, of fhe 
operafors Oa{Xa) and fhe values Vs{±q) faken by fhe so-called shiff function on fhe lefl/righl endpoinl of fhe Fermi 


/ 2n'\p(Vs(?)-Oj+/f)+p(Vs(-^)+;^) 


^ L > 


•e 


•Tl(-v,(-^), -k I OJs)yR{Vs{q)-Os,K-Os \ tu^) (0.2) 
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zone. A more detailed description of the correspondence and the constituents of the operator ffs can be found in 
the core of the paper. We do stress that the expectation value appearing in the rhs of (10.11) is trivial to estimates 
by means of the free fermion exchange relations, cf. (ID.5I) - (ID.6I) . We do stress that the correspondence (I0.1I) - (I0.2I) 
only provides one with the leading order to each oscillatory harmonics. This leading order only involves primary 
operators. However, should one be interested in sub-leading corrections then one would need to supplement the 
expansion (10.21) with additional terms that would involve descendants. Obtaining these terms is, in principle, 
possible within the method developed in the present paper, although would demand more computations. 

On top of providing the microscopic origin of the appearance of a c = 1 conformal field theory as an effective 
theory in the large-distance regime, our setting supports the interpretation of the universality hypothesis that has 
been advocated by one of the authors in 1341. Namely, within our approach, the leading power-law contribution 
associated with a given oscillating harmonic in the large-distance asymptotic expansion of a correlator stems 
from a saddle-point like contribution that is extracted from the form factor series expansion of the correlator. In 
the large-distance regime studied in the present paper, the saddle-point is located on the two ends of the Fermi 
zone. In the time-dependent case, one should as well incorporate all the possible critical points of the dressed 
momentum/energy combination xp - te, cf. 1261 . It is the singular structure, viz. local behaviour, of an operator’s 
form factor in the vicinity of the saddle-point that fixes fhe value of the critical exponents. Thus, two models 
belonging to the same universality class -in the sense that their appropriate correlators share the same critical 
exponents- have to share the same singularity structure of their form factors. Note that the regular part of the form 
factors can differ from one model to another in that that they solely impact the value of the amplitudes, viz. the non- 
universal part of the large-distance asymptotic expansion. Recall that the standard formulation of the universality 
hypothesis states that two models sharing the same overall symmetries share the same critical exponents. In this 
light, it would appear that the set of symmetries of a given model completely fixes the singularity structure of 
its form factors. However, it is not clear for us at the moment how one could derive the singular structure of a 
model’s form factors solely building on its symmetries. Nonetheless, it seems more satisfactory to us to think of 
universality in terms of classes of models sharing the same singular structure of their form factors. Indeed, this 
criterion is precise and can be checked, in certain cases, by explicit calculations, say in a perturbative regime. The 
matter is that there is, a priori, no criterion allowing one to say that one has identified all the symmetries of a 
model that are pertinent for fixing its universality class. It could well be that two models have apparently the same 
symmetry structure but that one has missed some important yet quite oblivious symmetry of one of the models 
that would imply that, in fact, the models belong to two distinct universality classes. It is worth, in this respect, to 
remind that two apparently very similar models, the 2D Ising and the the eight-vertex model in its lattice model 
formulation belong to fundamentally distinct universality classes. 

The paper is organised as follows. The present section is the introduction. In Section [T] we recall the free 
fermion based description of the space of states of the free boson model. We also present new formulae for the 
expectation values of exponents of specific current operators. In Section [2] we present the general framework 
-properties of the model’s spectrum and form factors- that allow us to derive the correspondence with the free 
boson model. Finally, in Section [3l we establish the main result of the paper, namely the set of formulae (10.11) - 
(10.21) that appeared in the introduction. The paper contains four appendices. In Appendix |3 we briefly review 
the special functions which arise in the course of our study. In Appendix |B] we compute some two-dimensional 
integrals which are of interest to the probelm. In Appendix O we prove Proposition II.II Finally, in Appendix 
IdI we show how one can recover, starting from the formalism developed in the present paper, the value of the 
multi-point restricted sums that were first introduced in I271I . 
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1 A free fermion description of the free boson model 

In the present section, we recall the free fermion based construction of the space of states for the free boson model. 
Our presentation basically follows the notations and conventions that can be found in the excellent review paper 
in . The various results found in this review originate from a long series of developments which started with the 
cornerstone work of Kyoto’s school (limbo, Miwa and Sato) in the late ’70’s on holonomic quantum fields l[2]l . 

1.1 Overall definitions and generalities 

1.1.1 The space of states 

We consider a set of fermionic operators {i/^n}«£Z and their * associates {i/^*)„ez which satisfy the anti-commutation 
relations 


{lAn.tAm) = {lAn><Am} ^ 0 and ^ Sn,m , (1-1) 

where is the Kronecker symbol. We also assume the existence of a vacuum vector 10 ) which satisfies fo fhe 
properties 

t}/„\Q) = 0 for n<0 and lAj^lO) = 0 for n>0. (1.2) 

The dual vacuum ( 01 fulfils fhe analogous properfies 

= 0 for n < 0 and = 0 for n>0. (1.3) 

Sfarfing from fhe vacuum (resp. fhe dual vacuum) one consfrucfs vecfors (resp. dual vecfors) fhrough a repefifive 
acfion of fhe fermion operafors. Such vecfors are paramefrised, in a nafural way, by fhe sefs 

Jn,-n, = [{Pa]7 \ [haf,^] (1-4) 

consisfing of fwo collecfions of ordered infegers 1 < pi < • • • < pn^ and 1 < /ji < • • • < If is convenienf, for a 
deeper physical insighf, fo fhink of fhe infegers {pa}”’’ as a labelling of fhe parficle-like excifafions and fo fhink of 
fhe infegers {/Jq}"'' as a label for fhe hole-like excifafions. To be more precise, fo each sef defined as above, 

one associafes fhe vector 


\J'np-,ni,) — 0 ) (1-5) 

and fhe dual vector 

{Jnp-,n„ I - < 0 ■ (1-6) 

The Hilberf space h of fhe model is fhen defined as fhe span of fhe vecfors infroduced above 

h = span||J„^;„^) wifh rip, tin e N and | < ... < e N*|. (1.7) 

Note fhaf, when fhe number of parficle and hole-like integers coincide, i.e. rip = nr, = n, one can identify fhe 
sef J'ripini, wifh a Young diagram. The one-fo-one map is obfained by inferprefing fhe infegers |{pa}” ; {/zq}”} as fhe 
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Frobenius coordinates of the Young diagranJll Such an identification is also possible when rip rih. The bijection 
is, however, more complicated. We refer the interested reader to the proof of Lemma [TTT] where one can find all 
of the necessary details for constructing such a bijection. Note that such a bijection was constructed for the first 
time, although indirectly, in Appendix A.3 of If24]| . 

The Hilbert space h also enjoys of another basis | Y; f) which makes the connection with a Young diagram Y 
explicit. This new basis singles out the collection of the so-call vacuum and dual vacuum states having a prescribed 
charge t. 


I • • - fAolO) f > 0 


and 


in 


^>0 

(0|(A-i • • • ^<0 


( 1 . 8 ) 


General states of this basis are build as equal in number particle-hole excitations over the vaccua | £ ), resp. their 
dual vaccua {i\. Let 


Y = ({aj” : {ySJ';) with 


( I < a\ < ■ ■ ■ < an 
1 < /?! < • • • < 


(1.9) 


be the Frobenius coordinates of the Young diagram Y. The second basis we mentioned above takes the form 

(Y;^| = 

In fact, one can even consider mixtures of the basis | Y; f’) and \J'np-,ni, ), namely the basis 

\3np-,nh'^£) = • • • <Api+r-iK) ■ (1-11) 

It is easy to convince oneself that 

Span{|Y;^-i-r) : Y Young diagram} ^ Span{| ) : sets np-nh = £ + r] 

= Spanjl : sets Up - nn = r\ . (1.12) 


1.1.2 The space of operators 

The fermionic operators lA/ and lA* can be thought of as building bricks allowing one to construct more general 
operators on the Hilbert space f)- The simplest class of operator is obtained as a linear combination in the fermions: 

m m 

where Vm, Wm are bounded sequences. The field and conjugated field operators 

T'(z) - and T^Cz) - (1-14) 

fsZ jeZ 

are an archetype of such operators . 

^we adopt the slightly unusual convention where the origin of the Frobenius coordinates is the diagonal so that these start from 1. 


7 



The second important class of operators is obtained by means of a bilinear pairing of the fermions realised by 
an infinite matrix A having bounded entries: 

(1-15) 

iJeZ 

Expressions of the type (11.131) or (11.151) may look a bit formal in particular due to convergence issues of the 
sums. However, all the sums of interest truncate as soon as one computes matrix elements taken between the 
fundamental system of basis vectors labelled by the sets J'n ■,nh ■ H is in this sense that all of the above expressions 
should be understood. 

The most fundamental example of an operator built through a bilinear pairing is the charge operator which 
takes the explicit form 


Jo = (1.16) 

k>0 k<0 

It is readily seen that the vector | J'n ) is associated with the eigenvalue np - n/, of the 
will sometimes say that the vector \J'np;ni, ) has charge np - «/,. 

In other words, the charge operator induces a grading of the Hilbert space f) into the 
having a fixed charge J 

h ^ with = spanjl ) : Up - nu = i]. (1.17) 

teZ 

The current operators provide one with another example of important operators that are bilinear in the fermions. 
To each half-infinite sequence t - {4}jt€N* one associates the current operators 

J±{t) ^ Y^tk-J±k with Jk = foi" ^^0- (1.18) 

k>\ jeZ 

It is readily checked that the components of the current operator satisfy the algebra 

[Jk,Je\ = k6k,-e . (1.19) 

These commutation relations allow one to interpret the currents Jk as bosonic operator modes (see e.g. lUll. 
Furthermore since, for k e N*, one has 


charge operator Jo- We 
direct sum of spaces hf 


7i|0) - 0 and (0|7_^ - 0 


( 1 . 20 ) 


it follows that one can interpret {Jk}kew ns bosonic creation operators while {/-jtlieN* ns the bosonic annihilation 
operators. It is on the basis of such an observation that the construction of the free boson field theory is made in 
the free fermion model we have introduced so far. 

It is also straightforward to check the commutation relations 

[Jk,4'i] = ^t-k , [Jk,^*e\ - -^l^uk ■ ( 1 . 21 ) 


The above commutation relations can be conveniently encoded on the level of the field and conjugated field 
operators. As a matter of fact, then, they take the form of the following exchange relations 


T(z). ^ 


and 


T*(z). e 


J±(t) 




( 1 . 22 ) 






where ^{t,z) = T,k>i Given a generic sequence t, this is the most compact form of the exchange relation 

since the sum defining ^ cannot be computed explicitly. However, for the specific choice t = t± where 


(G)* 


k 


and 



(1.23) 


fhe sum reduces fo fhe Taylor expansion of fhe logarifhm. If is fhe t - t+ choice fhaf plays a parficularly imporfanf 
role in our analysis in fhaf if is associafed wifh fhe consfrucfion of verfex operafors. We fherefore infroduce a 
special nofafion for fhe currenf operators associafed wifh fhe paramefers t±, namely 

- J±{t±) = J±k- (1-24) 

k>l 

In fhe case of fhe operafors ^±(v, a>), fhe exchange relations (11.221) particularise fo 

and = (l-(£)*^)-£^*ti,^).'F*(z) . (1.25) 

There is also anofher operator fhaf plays a role in our selling, fhe so-called shifl operator which we shall denote 
as e^. The shifl operafor maps hr onlo hr+i ■ In fncl, fhe operalor P arising in fhe exponenl is fhe conjugale operalor 
to 7o- Allhough fhe shifl operator has no simple expression in terms of the fermions, it takes a particularly simple 
form in the basis of h subordinate to Young diagrams, cf. (11.101) . Indeed, any integer power e''^ of the shift operator 
satisfies e''^| Y;^) = |Y;f?-i-r), viz- 


e'-^ - 2]|Y;^ + r)<Y;^| (1.26) 

Y,e 

where fhe sum runs over all inlegers £ and all Young diagramms Y( see |T| for more delails). 

We end Ihis sub-seclion by introducing the r-shifted bosonic vertex operators which play a central role in the 
correspondence: 

r(y, rim)- . ^rP _ ^^ 27 ) 


1.1.3 The Wick theorem 

In this sub-section we review the statement of Wick’s theorem in the case of an insertion of group-like elements. 
This name refers to a class of operators G satisfying to the so-called basic bilinear condition: 

2] (g ■ ^j) ® (g ■ ri) = 2 (<Ai • g) ® • g) . (1.28) 

;eZ ;eZ 

Group-like elements single out, among other things, because they allow one for a powerful generalisation of the 
Wick theorem which we shall recall at the end of this sub-section. 

The operators we have introduced so far provide several examples of group-like elements. For instance, it is 
easily checked that the operators eAd) in general, jn particular, and e'"^ all satisfy to the basic bilinear 

condition. It is likewise trivial to establish that, if G,G' are group-like elements, then so is their product G ■ G'. 
These facts ensure that that the r-shifted vertex operators Y'(y, r | a>) are group-like elements. 

We stress that group-like elements need not to be invertible. For instance, given any group element G and a 
collection oi,..., d„, resp. ro*,..., ro*, of linear combinations of the basic fermion operators tf/j., resp. i//*, viz- 

Vk ^ ^Vkm^m resp. wl = (1-29) 
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the operators 

-G and = G • ro* • • • ro* d„ • • • oi (1.30) 

are also a group-like. We shall not enter deeper into the detail here and refer the reader to fT] for more details. 
Group-like elements allow for the following generalisation of Wick’s theorem. 


Theorem 1.1 Let vj, resp. ro* , j = . ,nbe any operators linear in the fermions and G, G' group-like elements. 

Then the associated expectation value admits a determinant representation: 


<0|Gu)*---tn>„---DiG'|0) 

<0|G-G'|0) 


det„ 


r<0|Gu);t),G'|0)i 

<0|G-G'|0) 


(1.31) 


1.2 Expectation values of exponents of current operators 

In the present sub-section, we compute the expectation value of the operators between any two 

states I J'np,ni, ) I fTfu.n, ) of the basis. The corresponding result will provide a first brick towards the correspondence 
with the free boson model that will be built in Section [3] 

Prior to stating the result, we first need to introduce a few handy notations. Given two sets of integers 


Jnp,n,. - {{Paf{ \ and ; Ua}"'] 

one defines the set functions 


nt 




ni, (Yii^-h-ha-^v) n„ .Y\{pa + tb + V- 1 ) 

= f[- -n 


_ _ i i 

Y\{tb-ha-^v) a=\ Y\{Pa-kb-t- v) 

b=\ b=\ 


and 


(1.32) 


(1.33) 


< 3=1 ^ 




Pa 


a=l 


X 


hpt 


Up III, 

n {Pb - Pa) - n ih - ha) 

a>b a>b 

nh 

n Y\{Pa + hb- 1) 

a= \ b=\ 


(1.34) 


These functions appear as the building blocks of the so-called discrete form factor 




(jnp-,n,\Jn,-,n,\v,C0) - (-1)«P+”' • (-1) 


/sinfTTv] 


I . 


^{Jnp-n, I V,m) 


X I -V,m • m{jr,^.n,y,Jni^.n, I v) . (1.35) 


In (11.351) we made use of the hypergeometric-like notations for ratios of T-functions, cf. Appendix lAl 
It is the above discrete form factor that enters in the description of the expectation values of the current 
operators taken between the states \J'nk\n, ) and {J'np;ni, I- 
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Proposition 1.1 Let Snp\nh J'nf,n, be as given in (I1.32I) . Then form factor of the product of exponents of 
current operators admits the following representation 




'p-nh,nic-ni 


•^{fTnp,ni,ifTni;,n, I • 


(1.36) 


Note that the Kronecker symbol in (11.361) is there to emphasise that the form factor is zero unless np - n^ = 
nt - nt- This is due to the fact that the exponents of current operators preserve the charge of a state, viz. can only 
connect states belonging to the same charge sector. 

The above statement is well known in the case np = n/, = 0. Indeed, then, the representation (11.361) follows, 
for instance, from the Jacobi-Trudi determinant representation for Schur functions along with the the possibility 
to compute explicitly the Schur functions Sap{t±) associated with hook diagrams (a,/3) in Frobenius notations (c/. 
Appendix A of [IJ for more details). We do however stress that the general case given above is, to the best of 
our knowledge, new. In particular, the proof of the general case relies on new ideas not related to the theory of 
Schur functions. The proof, being slightly technical, is postponed to Appendix O We would also like to point out 
that Proposition 11.11 yields a new type of explicit representation for the skew Schur functions associated with the 
parameters t+. 

In order to identify the emergence of an effective field theory description in the large-distance spacing between 
the operators, one needs to determine the matrix elements of the general r-shifted vertex operators that have been 
defined in (11.271) . These will be given in Theorem ll.2l below and can be deduced from (11.361) . 


Theorem 1.2 The form factor of the r-shifted vertex operator Y{v,r\ co) reads 


{Jnp-,njny,r\aj)\Jn,-n,) = (- 1 )“ 


rlr+l) dfi —nh,nii—n,+r 


G 


(m) 


!^+r{nk-n,) \ l-V-r 


1 -y 


^{jnp-ni,',Jnf,n, I V, w) . (1.37) 


The symbol G appearing above stands for the hypergeometric-like notation for the ratio of two Barnes functions, 
see Appendix^ 


The proof of the theorem builds, in particular, on the form of the action of the r-shifted shift operator e''^ on 
the states I27n*;«, )■ We shall determine the form of this action in Lemma [TT] given below. Still, first of all, we 
need to introduce some notations relative to so-called r-translates , of the set 27n*;n,- 

Given a set fTn^-n,, we define q, 0 < q < r, tobe the unique integer such that 


1 < ti < • • • < < r < < • • • < . 

Then, the set {ta}\ ^ stands for the complement of the set {r -i- 1 - in [[ 
1 < 7i < • • • < Tr_q < r and {r -i- 1 - G)i U = 11 

The r-translate , of the set 27n*;n, corresponds to the set 


TA = [»'Si ■. 


; rj, namely 
rj . 


The sequence k'p^ appearing in (11.401) reads 

k' = ta for a = 1,... r - <7 and - ka -v r 

and has a total of = n^ + r - q elements while the sequence t' is defined 
t'a = ta+q-r for a ^ I,.. .n, - q 


for 

as 


a - . .nk 


(1.38) 

(1.39) 

(1.40) 

(1.41) 

(1.42) 
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and has a total of n' = rit - q elements. 

Proof — (of Theorem M .2\l Let r > 0. Then the representation (11,371) follows from the form taken by the action of 
the shift operator (11.451) on the states \ffnk-,n, ) and from the identity 




- G 


( • ^(jn,-nP,Jn,-n, \v-r,aj). (1.43) 


There q is as defined in (11.381) while (-l)‘^4lfal) j-gfers to the signature of the permutation introduced in (11.461) . 

The formula (11.371) holds, in fact, for r < 0 as well. To establish this fact, observe that the whole construction 
is invariant under the transformation (tA/t, ^ It is easy to see by using the properties of the 

states and operators involved that, under this transformation, one has 

\'ny, r I Oj)\J„,-n, ) (- 1)”^-”'^'^^^"--"''^ I^C-V, -r | aj)\Jn„n, ) ■ (1-44) 

After some algebra, one can convince oneself that the rhs of (11.371) satisfies to the same transformations. 


Lemma 1.1 For any r > 0, it holds 

e'-'^ • \Jn,,n, ) - • (-1)-4IU). ) _ 

Above, CriUa]) corresponds to the permutation 

CTriita]) : I,-- - ,r 1-^ + I,-- - , r-ti + 1, tr-q,--- , ti 

(_l)o-r(hfll) Stands for its signature, while q and [ta) are as defined in (ll.38l) - (ll.39l ) 


(1.45) 


(1.46) 


Proof — 

We shall prove the lemma in the case when n^ - nt = i > 0. The 7’ < 0 case follows from a similar reasoning 
so that its proof is left as an exercise. Let s € [[ 0; 7 ]] be the unique integer such that 

1 < ki < ••• < < 7 < ks+\ < ■■■ < k„i^ . (1-47) 

We then introduce the complementary sequence {£ + \ - kaYf^ to {ka}\ in H 1; 7]]. The sequence {kaYf^ is such 
that 

1 <Ii < ••• <I^_, < 7 and {kJ^U{7+l-I«lf-^ - El;73 
The sequence {a^ is then defined as 

a a - ks+a - 7 for a = I, s 

and the sequence {/3a as 

/3a - ka for a = I,... ,f - s and /3e-s+a - ta + i for a = 1, ...,«? . (1.50) 


(1.48) 

(1.49) 
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Note that there is, in total, n^- s labels for the /3’s due to the condition tit - Finally, for convenience, we 

introduce the sequence ^ which we define = ( + \ - kfs+x-a^ so that one has the ordering 

1 < < ••• < < t. ( 1 . 51 ) 

Defining o^t{{ka}) as fhe permufafion 

ae{[ka}) ■■ ^ ks,... ,kuk\,... ,k'i_^ . ( 1 . 52 ) 

one has 

( 1 . 53 ) 

We denote by Y fhe Young diagram wifh Frobenius coordinates |{o'a)j*'~'' ; {ySa)”*~'*|. The vector | Y; ) can fhen 
be recasf as 


|Y;0 - (-1) 


o-e{\ka\)+ 


2 


-i • • • *^-'1 ■ ■ ■ • • • (Atri-rlA-tJ-l • • • 0 ) 

^ ( _ ifM . (_ 1 )^ . ( _ lf,+n,W-s) . ^ 


We now recasf fhe vector | Y; C + r) where fhe parfifion Y is defined in ferms of exacfly fhe same Frobenius 
coordinafes. Permufing fhe fermionic operators labelled by fhe indicies {r + I,... ,r + £} and {1,...,r) we are led 
to 


r(r-l) 


— m l^ , <(<-!) 

|Y;^ + r) - ( - . ( - 1)-^ 


+0'r(l'fll) 


/* /* /* /* 
• ■ *Ar-n ' ' ' 


X >kr+k„i^-l ■ ■ ■ <Ar+A:.,+i-l ' <Ar+t:,-l ' ' ' <Ar+t:i-l ' fAr+iJ-l ' ' ' ' ' ' <A7i-ll 0 ) 

_ ^ ^ ^ _ ^•^(n,+nk)(e-s)+q{n,+nk-g) , ^(r) 

The acfion (11.451) follows sfraighlforwardly from fhe above and from fhe represenfafion (11.261) of fhe shiff operator. 


I4«>- 0-55) 


2 The model 

In fhe presenf section, we describe fhe selling in which we assume fhe one-dimensional quanfum Hamillonian to 
til in. More precisely, we describe fhe assumpfions we make on ils speclrum, ils slates and on cerlain general 
fealures salisfied by fhe form factors of fhe model’s local operators. Such a description is sfandard for quanfum 
infegrable systems in one-dimension, see e.g. |8 ] in whal concerns fhe speclrum par! and fhe papers ll23l |25l in 
whal concerns fhe form factors. The validify of fhe piclure has been also esfablished, on a perfurbalive level, in 
i45] l. We Ihus believe lhal Ihis setting is common, in particular, lo all models belonging lo fhe universality class of 
fhe Lullinger model. 

2.1 General hypothesis on the spectrum 

We consider a one-dimensional quanfum Hamiltonian Y/ representing a physical syslem in finite volume L. The 
volume L represenfs, for inslance, fhe number of siles in fhe case of a lattice model or fhe overall volume of fhe 
occupied space in fhe case of a model already formulafed in fhe continuum. 
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We shall assume that the eigenstates of the Hamiltonian can be organised into sectors with a fixed pseudo¬ 
particle number N'. In practical situations, the integer N' may be related to the total longitudinal spin of a state 
(in the case of spin chains with total longitudinal spin conservation) or may simply correspond to the number of 
bare particles building up the many body eigenstate of (in the case of models enjoying a conservation of the 
total number of bare particles). 

In such a setting, the ground state of the model | '1'^;^? ) is located within the sector with N pseudo-particles. 
This number does depend on L and is such that, in the thermodynamic limit L —> -i-oo, one has limi^+oo {NIL) = 
D > 0. Furthermore, taking the thermodynamic limit restricts the space of states to the sector corresponding to 
excitations having n finite, when L —> -i-cx), energy relatively to the ground state. We assume that the eigenstates 
having this property are located in sectors with N' bare particles where N' is such that the difference s - N' - N 
remains finite in the thermodynamic limit. Having in mind the CFT-based interpretation of the integer s, just as 
the spin-chain setting, we shall refer to s as the spin of the excited state. 

We further assume that the excited states are only built up from particle-hole excitations. This means that we 
assume that it is possible to label the eigenstates, within each sector built up from N + s quasi-particles, by a set 
of integers 


I 


(^) 

n 




( 2 . 1 ) 


containing two collections of integers which label the so-called particle and hole [h^a^Y{ excitations. In this 

parametrisation, the integer n may run through 0,1,... -i- 5 while : 


< • • • < and < • • • < with 




. ( 2 . 2 ) 


The precise values of the integers defining the range of the p^’s vary from one model to another. Typically 
for models having no upper bound on their energy, one has = -i-oo while for model having an upper bound, 
are both finite but such that -N, a = 1,2, both go to -i-oo sufficiently fast with L. According to this 
setting, we shall denote the eigenstates of the model as 1). 

The vectors | ) provide one with the so-called microscopic description of the model, namely a complete 

parametrisation of the space of states in terms of discrete integers. However, it is the macroscopic description that 
is pertinent for describing the thermodynamic limit of the observables in the model. This macroscopic description 
arises by means of the so-called counting function associated with each given excited state ). More 
precisely, the particle, resp. the hole-like excitations in a given eigenstates are described directly by a set of 
rapidities {p^^J}”, resp. These rapidities are defined as the unique solutions to 



and 


j. /—(i)\ _ 


(2.3) 


Since, as mentioned, the counting function does depend, a priori on the set of integers labelling the excited state 
I ), the system of equations (12.31) is, in fact, extremely involved. 

In the following, we shall build on the assumption that any counting function admits, in the L ^ +oo limit, 
the asymptotic expansions 


^p.^){co) ^ ^(OJ) + 7^-i(m) 


L " 




(2.4) 


The asymptotic expansion of the counting functions involves three "macroscopic" functions. 
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The function ^ is the asymptotic counting function. It is the same for all excited states and we assume that 
it is strictly increasing. This function defines a set of "macroscopic" rapidities {/lalagz 


= - . 


(2.5) 


These macroscopic rapidities provide one with the leading order in L approximation of the rapidities 


and 


—(i) 


and 


-~<.v) 




( 2 . 6 ) 


• The function F^(s) stands for the shift function (of the given excited state in respect to the model’s ground 
State). It is a function of the spin .y and of the set macroscopic rapidities 


_ 

/\n — 




(2.7) 


This function measures the small 0(L ^) drift in the position of a rapidity in the Fermi sea under the effect 
of interactions. 


• Finally, the function represents the 1/L corrections to the ground state’s counting functions. It is this 
quantity that drives the non-trivial part of the first sub-leading corrections to the ground state’s energy. We 
stress that our way of decomposing the 1/L corrections to the counting function is such that the ground 
state’s shift function vanishes, i.e. = 0. 

In the large-L limit and within such a setting, the rapidities for the ground state form a dense distribution 
on {-q',q\ - the so-called Fermi zone of the model - with density . The endpoints +q are called the Fermi 
boundaries. The particle and hole excitations in the modefl carry a dressed momentum p and a dressed energy s. 
The dressed energy and momentum are smooth and satisfy to the general properties 

l[r> 0 e]_^;^[<0 and £r\[_^;^]>0. (2.8) 

The relative momentum and energy of an excited state are expressed in terms of the dressed momentum and energy 
as 


A£(ji^^) 


n 

13=1 

o 

+ 

(2.9) 


^ - nif) 

n 

a=\ 

o 

+ 

(2.10) 


Above, = {0; 0) refers to the set of integers which parametrises the ground state of the model and !P(J^^^) 

and are respectively the momentum and energy of the state parametrised by the set of integers 1^^^. We 

remind that the superscript {s) refers to the fact that the excitation labelled by takes place above the lowest 
lying energy level in the spin s sector. 

The excitation momentum and energy are probably the best examples allowing one to discuss a loss of infor¬ 
mation due to parametrising the particles and holes in terms of their "macroscopic" momenta. Indeed, on the level 

^Typically, for quantum integrable models, these functions are given as solutions to linear integral equations da. Their description, 
for more complex, non-integrable models, is in principle much more complicated. 
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of (I2.9I) - (I2.10I) . there arises a huge degeneracy in what concerns the 0(1) finite parts of the excitation energy and 
momentum in that numerous choices of the set will lead to the same value of the finite part. Indeed, the very 
definition (12.51) of the macroscopic rapidities ensures that 

^k-^k+i = O(^) • (2-11) 

Thus, a direct computation shows that for any fixed collection {ka, of L-independent integers, the two sets of 
integers 

In'^ = [{Pa}1\ mi] and = [[Pa+ka}1\ [ha + ta\l\ (2.12) 

will generate the same finite part of A£ or A!P, namely 

A£(ji^)) - A£(J(^)) - 0(i) and = 0(i) . (2.13) 

Thus, the so-called microscopic variables allow one to distinguish between all of the states of the Hilbert 
space whereas the macroscopic variables, at least in what concerns the leading in L order, naturally give rise to 
a huge indeterminacy in identifying a given state in the model - many microscopically distinct states give rise to 
the same macroscopic variables-. Hence, in the large L limit, there is a loss of information when passing from the 
microscopic to the macroscopic description. 


2.2 The operators and their form factors 


We assume that the model is endowed with a collection of local operators O^- These operators are best charac¬ 
terised in terms of their form factor, namely their expectation values taken between two excited states. We shall 
assume that the operator Or only connects those eigenstates which differ by Or in their spin, namely 


I x) = 


- e 


ixmim 




<jL''^ia(0)|J^'O ^ 0 only if 5-/ 


r{^) 


= Or 


(2.14) 


Scalar observables introduced in the last subsection were parametrised, in the large-L limit, solely by the 
macroscopic set of rapidities subordinate to the set of multi-indices This is no longer the case for form 
factors as we demonstrated in our previous work |[25]| . Within our setting, the latter are parametrised, in the large-L 
limit, not only by the sets of macroscopic rapidities but also by the sets of discrete integers 

Namely, for properly normalised states | ) and their duals ('P(J^^) |, the form factors take the form 


ToXim ; 10 ) = 


(A. 

■m 




,(i+Or) 




A). 

•m f 


r:, 


(i+Or)| j-is) . j-is+Or) 
m rt 


’)■( 


1 -F 




(2.15) 


In this decomposition 

• is called the smooth part: it solely depends on the macroscopic rapidities ,R^Xi ■ Furthermore, 

this dependence is smooth what ensures that a change in the value of the integers in the spirit of (12.121) will 
only affect the value of the smooth part in the 1 /L corrections. 


• The smooth part is a set function so that it is invariant under permutation of the particle or hole rapidities 
meaning that 


with 

I’k "l 

The same type of reduction holds as well for particle-hole rapidities 
spin s + Or sector. 


«!? = ItMkpU ■ ■ 

associated with the excited states in the 
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is called the discrete part because it depends explicitly on both types of parametrisations: the macro¬ 
scopic rapidities and the sets of integers labelling the excited states We do stress 

that this last dependence is explicit in the sense that it goes beyond a sole dependence of the integers through 
the parametrisation by the macroscopic rapidities. The main effect of such a dependence is that a change in 
the value of the integers in the spirit of (12.121) will result in a significant change (of the order of 0(1)) in the 
value of The discrete part thus keeps track of the microscopic details of the different excited states. 

The smooth part represents, in fact, a non-universal part of the model’s form factors. Its explicit expressions not 
only depends on the operator Or but also varies strongly from one model to another, see Il25l |33l for concrete 
examples of the form taken by Sq^ for different quantum integrable models. The part is, however, entirely 
universal within the present setting of the description of the model’s spectrum. It solely depends on the values of 
the spin Or of the operator Or- Its general explicit expression plays no role in our analysis, in the sense that we 
shall only need the expression for specific excited states, namely the ones belonging to the so-called /’^-critical 
classes that we define bellow. 

We refer the interested reader to ||2^ for a thorougher discussion relative to the origin of the discrete part 
within the framework of Bethe Ansatz solvable models. 


2.3 The critical 4 class 

Observe that, in virtue of (12.81) . one has s{q) = 0. As a consequence, by (12.91 ). there exists a possibility to realise 
zero energy excitations in the thermodynamic limit L —> -i-oo by forming holes and particles whose rapidities scale 
down to ±q when L —> -)-oo. This singles out a class of excited states of the model which we shall refer to as 
critical states. An excited state is said to belong to the critical class if , in the L —> -)-oo limit, all the macroscopic 
rapidities of the particles and holes associated with this state "collapse" on the Fermi boundary: 


pts) - ±q and - ±q . (2.16) 

Pa '^a 

There, the ± sign depends on whether the particle or hole’s rapidity collapses on the right or the left Fermi 
boundary. 

A set of integers is said to parametrise a critical excited state if the associated particle-hole integers [p^a^Yl 
and can be represented as 


- iAi ^ A9 t”p;+ I I (1 _ „6) 


[p'aXi = [N+s+^ [I - 

where the integers Pa^±,h^ali ^ "small" compared to L, i.e. 


lim - lim -^ = 0, 


L—»+oo L^+oo 

and the integers satisfy to the constraint 


(^0 , ( s ) (‘^) , (‘^) 

n J. + n j = n) ^ n) = n . 
P^+ p- /?;+ h\- 


(s) (s) 

+ U (2.17) 

(2.18) 

(2.19) 


Within this setting, one can readily check that the critical excited state described above will have particles, 

resp. holes, on the right/left end of the Fermi zone [-q ; g] associated with the spin s sector. 

In fact, one can distinguish between various critical states by organising them into so-called ^.^-critical classes. 
This classification takes its origin in the fact that all such states have a vanishing excitation energy (up to 0(1/L) 
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corrections) but can be gathered into classes depending on the value of their macroscopic momenta 2€sPf, where 
Pf = p{q) is the so-called Fermi momentum and 


(.v; 

p; 


,(h 

'h 





( 2 . 20 ) 


Note that the subscript 5 in allows one to localise the given critical class in a specific spin 5 sector. 

The finite part of the relative excitation momentum associated with an excited state belonging to the £s critical 
class only depends on £s, while the first dependence on the specific representative of the class arises on the level 
of the corrections. More precisely, one has: 


APiln^) = IfsPF + - 1) + “ ycr I ^{Pa - - 1) + 


«/t;+ 


a=l 


a=l 


In 


nir- 


2n ( _ ^ 
+ —[a I 


a=l 

{£s + 1 ) 


< 3=1 

+ . (4-1) 




+ ... ( 2 . 21 ) 


where we have set 

^ 1 p'i±q) 

2n^'{±q) ■ 

We insist that all the terms included in the dots either 


( 2 . 22 ) 


• are of the order of 0(1/L) but do not depend on the integers pa-± and ha-+ nor on np/h-+ (they can, nonetheless 
depend on is)', 

• depend on these integers but are of the order of 0{\/L?). 

Note that, in principle, the terms present in the second line of (12.211) could have been simply included in the .... 
We chose to write them down for normalisation purposes. Namely, the formula (12.211) in its present form allows 
one for a more straightforward correspondence with the free boson model. 

The partitioning of the particle and hole integers of a given state belonging to the is critical class into two 
collections of integers of ± types (12.171) suggests to parametrise the excited states within the is critical class 
directly in terms of the two subsets . U . where 


rris) 
•J n\m 



Hence, from now on, for any state belonging to a critical is class, we shall identify the sets 


I 


(h 

n 


and 


j(s) 


(.9) 


(s) 

*;+ 


u 


(•?) 

n 

rtp- 


\n 


(S) • 

In- 


(2.23) 


(2.24) 


We do stress that, due to (12.171) . the value of the spin s does play a role in the correspondence between and 
() U . . We also draw the reader’s attention to the fact that the value of is is encoded in the very 

4'- 

notation U , see (|2.20D. 

444;+ 4;-4;i 
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2.4 Large-L expansion of form factors connecting critical states 


We are finally able to discuss the expression for the form factors of local operators O, taken between two eigen¬ 
states belonging to critical classes. It is precisely these form factors that are responsible for the emergence of a 
effective field theory description at large distances of separation between the operators. Let 


j(i) _ q- 
m — U f 




u j; 




and J 


(i+Or) 






be two sets of integer^ parametrising excited states belonging to the 

and ^in = np-+ - nh-+ = nu-- - Up-. 


(2.25) 


(2.26) 


f out ^ mp-^ - mh-+ = nih-- - nip-- 

critical classes. 

Within the framework we impose on the model, the form factors of local operators taken between two excited 
states belonging to the critical classes introduced above take the form: 


I Xr) = • Te^^,-e,SOr) 


/ '*'Lut An)'rp(t^,-+^out An) 

T/ 


X ^ 


[ 27 inp;+;m*;+ J 27” np-+-,ni,-+ I ] ' =^[27”i 27 np --,ni, - I ] 

7 -1 / 

_i_ / in '■ p ^oui 

(m^) ■" 2 ^ 

The constituents of the above formula are parametrised by the values 

= Vrici) - Or and V~ - Vr{-q) 

that the relative shift function between the 7in, 7out critical states 
Vr{X) = F,(A) - F.+. td) 


(m;) 2 2 


• (1 + 0{^fj . (2.27) 


(2.28) 


(2.29) 


takes on the right/left endpoints of the Fermi zone, up to subtracting the level Or of the operator Op in the case of 
the right endpoint. The also depend on 

-2ma- q- 


Ojt = e2‘™"T 


and 


o)^ - e 


(2.30) 


representing the exponent of the individual momentum brought by a particle excitation on the left or right Fermi 
boundary. In the above large-L asymptotics, the quantity represents the properly normalised finite and 

non-universal {i.e. model and operator dependent) part of the large-L behaviour of the form factor of the operator 
Or taken between fundamental representatives of the Cm and 7out critical classes. More precisely, it is defined as 


15;) 


J, 4n)"^'p(4 "'"^out An) . , 

\ / 


<4'^ io,(0)ix4'’''A 


(2.31) 


in which the sets of integers IsP and parametrising the excited states correspond to the fundamental 

CQut ^in 

representatives of the 7out and Cm critical classes. Namely, is the set of particle-hole integers living on the 
Fermi boundary in the spin s sector and such that 


r{s) _ 
— 


j{p4 = a\\ ; {0}) U {(01 ; (4- - ci]\] if 7 > 0 
j{0) ; {/j4 = a)~^] U \{ppi = a}-^ ; {0}) if 7 < 0 


(2.32) 


Tn this writing we chose not to omit the spin sector label so as to lighten the formulae. Since the context of the setting is clear, this 
ought not lead to any confusion. 
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Note that this form factor solely depends on the difference {out - An- This issues from the fact that, a priori, 
in the thermodynamic limit, any fundamental representative can be chose to be the reference ground state. Yet, 
relatively to the fundamental representative of the ^out-class, the other excited state corresponds to the fundamental 
representative of the {out - An class. 

The power of the volume L arising in (12.271) and (12.311) involves the right pCv^) and left p{v~) scaling dimen¬ 
sions whose generic expression reads 

p(y) = - . (2.33) 

The factors ^ have been defined in (11.351) and contain all the "microscopic" contributions issuing from ex¬ 
citations on the right or left Fermi boundary. These local microscopic form factors depend on the value taken 
on the right or left Fermi boundary by the relative shift function Vr associated with the critical excited states of 
interest and on the position of the operator. These local microscopic form factors also depend on the sets of in¬ 
tegers ±;m;, ± and 27np ±;n;, ± parametrising the excitations on the boundary of the Fermi zone to which they are 
associated. Those microscopic form factors are as defined in (11.351) . 

Finally, the normalisation constant v~) is chosen in such a way that it cancels out the L-independent 

contributions of the right and left critical form factors when focusing on the fundamental representative of the f’out 
and An critical classes. Its explicit form can be computed in terms of the Barnes G-function O and reads: 


C 


(^out All)/ 


(v^v;) = G. 


1 -H y^, 1 - y^ 

1 -)- y,. -H {out ~ An^ 1 ~ y^ ~ f^out + An 
Above we have used the hypergeometric-like notations for ratios. These are explained in (IA.5I) . 


(2.34) 


2.5 General remarks on the scope of applicability of the model 


The structure of the spectrum and form factors that we discussed throughout this section already appeared in the 
paper ||27]| . 

The formula for the form factors of local operators taken between excited states belonging to critical classes 
can be proven within the framework of the algebraic Bethe Ansatz for various quantum integrable models on the 
basis of determinant representations for their form factors ll28l 1^ 1431 . The corresponding calculations 

are a straightforward generalisation of the method developed in ll2^ |25l HQ. However, we do trust that the 
decomposition (12.27b is, in fact, universal. More precisely, the properly normalised form factor 
is definitely model dependent and thus can only be obtained on the basis of exact computations. Its explicit 
expression is available for many quantum integrable models. However, we do trust that the local microscopic 
pre-factors and the leading power-law behaviour in L is universal: namely that they take the same form for models 
belonging to the Luttinger liquid universality class. This has been confirmed by perfurbafive calculafions around 
a free model ||431 . 

The affenfive reader mighf observe a slighf difference befween fhe expression for fhe form facfors faken be- 
fween crifical excifed sfafes as given in Secfion IT4l and Section 2.4 of paper lITTl . More precisely, fhe formulae 
differ by 


i) an overall sign 


(- 1 ) 




(- 1 )^ 


■(«*;+-!) . 


(2.35) 


ii) in llT7ll . fhe expression for fhe normalisafion consfanf ^“^(y^, y^) look fhe form 


C 


(^outiAn)/ 

original 


(y+,y;) - G 


1 -)- y„ 1 


, 1 + Ain 


y^, 1 


1-A 


m + Vy,\ + 


■> I ^out + A 


■ An + y, 

y^, 1 -)- A 


r 

out 


Ain + Vy 


(2.36) 
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iii) the presence, in ll27]l . of the factor 


^ extra Vr’ 


sin[7ry;!']\4ut 
sin[7rv7] ' 


(2.37) 


in the definition of the local microscopic form factors as compared to the definition (12.271) used in the present 
paper. 


The difference of sign pre-fractors is irrelevant in that the sign can be readily absorbed into the normalisation of 
the eigenstates. The present choice of the sign was more natural here since it makes the comparison with the free 
boson model more direct. Finally, in what concerns the normalisation constants, one can check, on the basis of 
identity (IA.4I) . that one has 






_ /^(^out An)/ 


(2.38) 


3 The resurgence of an effective free boson field theory in the large-distance 
regime. 

3.1 An effective free boson field theory model 

The effective Hilbert space is defined as the tensor product of two copies hi and hi? of the free boson Hilbert space 
f) introduced previously: 


f)eff = f)L ® t)i? • 


(3.1) 


The first space (resp. second) arising in the tensor product will be called left (resp. right) space and will be 
associated to modelling what happens on the left (resp. right) Fermi boundaries. 

To a local operator Or{Xr), we associate the below operator on heff 

ffriOJr) = 2]r.(a)- 




In this formula, 

• '^lir{v, k', a>) stands for the operator acting non-trivially on the L/R copy of the original Hilbert space as the 
vertex operator y (y, x; co) defined in (11.271) . viz- 

Tz,(y, A-lm) = y{v, K \ (jS) ® id and "^(y,/c-|a») = id®y{v,K\(jo). (3.3) 


• p(y) are the scaling dimensions introduced in (12.331) : 

• yK{Or) is the properly normalised in the volume thermodynamic amplitude defined in (12.311) and taken 
between excited states satisfying €oa\. - An = k- 

• my is a phase factor that reads 


g±2i;r^a* 


(3.4) 
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It follows readily from the discussions carried so far that given two sets 


jis) _ q- 
m — t 




u j; 




and J 


(■«+Oi) 


^ J, 


«p;+;«/!;+ ^ 3'np --,nk - 


(3.5) 


parametrising critical excited states in the physical Hilbert space hphys, one has an equality, up to 0(lnL/L) 
corrections between matrix elements: 


{l ® I .5' + Or',Sfnp-+-,nt+ )} ' )) 


where 


>i = >np-+ - mh-+ - np-+ + nh-+ . (3.7) 

Note that, in order to establish (13.61) . one should use that 

{ j7mp;mA I “ |o and \ SFmp',mi, ^ ^ ^ ~ 0 ) (3.8) 

as readily inferred from Lemma [TT] 

Note also that (I3.5I) - (I3.6I) gives rise to a correspondence between states belonging to ^-classes in the original 
Hilbert space and general states in the effective Hilbert space: 

\Im) = U; Jnp,--ni,-,- ) ® U + O/, Jnp.y,n,,-p ) (3.9) 


3.2 The effective truncation of I)phys and the correspondence 

Consider the r-point ground-to-ground state expectation value {0\{x\)- ■ -Orixy))^ . One way to express this 
quantity is to write its form factor expansion, viz- insert the closure relation between each operator. According to 
our hypothesis on the space of states, this recast the correlator as 


r r-\ 

/ t{?s) 1 5=1 5=1 

Above we have introduced the sequence of spin gs - Zp=i Op and the summation runs over all choices of 
and sets with 5 =l,...,r-l. Note that we have adpoted the convention 

f„ = ?,. = 0 and 7“ = J<g, = ||»| ; |»|) . (3.11) 

The above sum is a highly oscillatory sum in the large-distance regime Ix^+i - Xs\pf » 1- Therefore, it shoud 
localise, in this regime, either around the saddle-points of the relative excitation momentum A!P(J^^^) defined in 
(12.101) or around fhe boundaries of fhe summation region. In the present setting, due to the hypothesis (12.81) . there 
are no saddle points, so that the leading in |x. 5 +i - Xs\pf » 1 contribution to the sums will issue from a vicinity 
of the boundaries of summation. This corresponds precisely to the case of the states belonging to critical classes. 
Thus, this reasoning leads to 


(1 s=\ ’ A;+ 




(fj-l) 

”p;- 


. rriFs) 




'(V.d 


,(?s) 

h-- 


(3.12) 


(f^) j 
'h-,± 
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Note that the - symbol refers to an equality up to subleading corrections to each oscillating factor 
with s = 1,..., r - 1 and is ^ Z. Hence invoking the identity (13.61) . using that 


Y]{-^r -1 

a=l 


with Xs 


m 


p;+ 


h; 


p; 


+ n 


(V.i) 

/i;+ 


(3.13) 


and ignoring the effect of the 0(lnL/L) corrections that should vanish in the thermodynamic limit, we re-absorb 
the sums over the intermediate states of I)eff^ hence obtaining 


(Oi(Xi) • • • - (^l(mi) • • • . (3.14) 

Thus the claimed correspondence holds. Above, we have added the subscripts hphys and heff so as to insist in which 
Hilbert space the expectation values are computed. 


3.3 Some examples 

3.3.1 The XXZ spin-1/2 chain 

We shall now illustrate the general framework described above on the example of the XXZ spin-1/2 chain embed¬ 
ded in an external magnetic field h > 0. This model corresponds fo fhe Hamiltonian 


l^xxz 


k=l 




V y 


k+\ 


+ A(o-^cr; 


k'^ k+i 



L 

with ^ cr| . 

k=\ 


(3.15) 


Above are the spin operators (Pauli matrices) acting on the k-th site of the chain, h is an external magnetic 
field and fhe model is subjecf fo periodic boundary conditions. If is well known fhaf 'Hxxz posses differenl phases 
depending on fhe value of fhe anisofropy paramefer A. When -1 < A < 1, ^xxz has a massless specfrum in fhe 
L - 1-00 limif. Below, we shall only consider Ibis regime and adopf fhe parametrisafion A = cos If is known 
fhaf, in ifs massless phase af /j > 0, fhe excifafions in fhe XXZ-chain can be eifher of bound sfafe nafure (so-called 
siring solutions) or be of particle-hole type. The parficle-hole specfrum enjoys of fhe slruclure described in fhe 
earlier pari of Ibis section. The bound slates, however, are described differenlly. None-lhe-less, one can argue 
fhaf, in fhe long-disfance regime, fhey only produce correcfions fhaf are exponentially small in fhe dislance; Ihey 
can fhus be disregarded when sludying fhe correspondence wifh fhe free boson model in fhaf if solely involve fhe 
algebraically decaying terms. 

Owing to [1-lxxz,S^] = 0, the role of quasi-particles is played by the number of down spins building up a 
given eigenstate. The local operators of the model consist in products of elementary matrices 


Okinik) 


mt+dt 


(3.16) 


The length dk of the string of elementary matrices is assumed fixed, viz. mk and L independenl. There is fhe 
elemenlary malrix acting non-frivially on fhe p* sife of fhe chain. The operafor defined by (13.161) carries fhe 
spin 

Ok = #(a : (e^;«,ty.J - (2,1)) - #(a : = (1,2)). (3.17) 

One can show using Befhe Ansafz mefhods fhaf fhe shiff functions - in fhe sense of (12.41) - fake fhe form 


F^(.,)(T) = s{^{A,q) - -t ^[(/.(4,//pJ - 0(T,p/,j], 

< 2=1 


(3.18) 
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in which 0 is the dressed phase, Z the dressed charge and the rapidities {iJ.a}aez are the unique solutions to 


a p(co) D 

^(Pa) = 7 With ^(w) = + — . (3.19) 

L 2n 2 

The function p arising in the expression of the asymptotic counting function is the so-called dressed momen¬ 
tum and D = limAr_»+oo(A^/T), N being the number of down spins in the model’s ground state. The function p 
corresponds to the unique odd solution to the linear integro-differential equation: 


<1 

pW + JK'I-P)- p'(p) ■ ^ 

-q 


/sinh(i^/2 -i- /l)\ 
Isinh(i^ 12- A)) 


(3.20) 


The functions Z and 0 solve the Lieb integral equations 


Z(/l) Jd'(A-p)Z(p)~ ^ 1 (p(A,v) + J0'(A-p)<p(p,v)~ ^ ^^’^( sinh(i7-d-ry) ) ' 

-q -q 

The above integral equation depend on the endpoint q of the model’s Fermi zone, which is fixed once that D is 
given. 

Finally, the relative shift function associated to critical excited states belonging to the and is-i classes 
and differing in Og in fheir pseudo-parficle numbers, fakes fhe form 

v.(d) = - 4>(A,q)) + (4-1 - ^.)(Z(T) - l) . (3.22) 

One can fhen show using cerfain idenfifies safisfied by Z(q) and (p{q, q) fhaf 


Vsiq) + K-Os 


= K^{q) - 


Os 

2Z{q) 


and 


Vs{-q) + K = KZ{q) - 


Os 

2Z{q) 


(3.23) 


We do sfress fhaf numerous laffice operafors Ok{nik) will carry fhe same spin Ok and fhus fhe same scaling dimen¬ 
sions since fhese are paramefrised by Vs{q) + k - Os and v^i-q) + k. 

Finally, fhe behaviour of fhe form factors of fhe model follows exacfly fhe form provided in Secfion|2] We do 
nol reproduce fhe corresponding formulae here buf refer fhe interested reader to fhe aforementioned liferafure. 


3.3.2 The Luttinger structure recovered 

One can reproduce the structure of the Luttinger liquid-like critical exponents by making additional assumptions 
on certain ingredients of our general framework. We first assume that the shift function F^(s){A) satisfies the 
reflection property: 


F^(s){-A) - s-F^(s){A) 


with = \{-p^(.)}\ ; {-/!/,w)^} - 


(3.24) 


This property can be argued, on the heuristic level, as follows. The physics should be invariant in respect to 
reflections of the Fermi surface. When taking the reflected of an excited state belonging to the N + s quasi-particle 
excitation sector above the ground state, one should reflect all the rapidities of the particles and holes but also take 
into account that the s additional quasi-particles that are used to construct this excited state which are initially 
located on the right boundary are now located on the left one. This induces the shift by -s. These properties can 
be directly verified on the example of the XXZ chain introduced earlier on. 
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Now, we shall make the additional assumption of a linear response, namely assume that if the macroscopic 

(3.25) 


variables ^ can be partitioned as 


(n{s+s') _ (u{s) I I 

R-n+n' - U 


then one also has 




(3.26) 


It is then enough to observe that the set of macroscopic rapidities of an i class, € >0 can be decomposed as 

_ f i 

Rf = [{q ]\; {-q]\ u) - U ^ U l<®* ^ ' ^1 


~^umkp 


='R« 


so as to get 

F^«(d) - + ^^«sp„(^) (3-28) 

Then the shift function arising in our estimates is recast as 


(3.29) 


where we have set k - Com - ^m- We do stress that all the above properties are verifies in quantum integrable 
models, the XXZ chain in particular. 

Our relations lead to 


Vr{-q) = v-j/cX-j) with vK - and j - F^^^^i-q) . (3.30) 

Then, the reflection property implies that Vriq) = Or - Vr{-q)- All-in-all, this reproduces the scaling dimensions 
of the Luttinger liquid model. 


Conclusion 

In the present paper we have provided a first principle-based derivation of the emergence, in the large-separation 
regime between the operators, of an effective description of the model’s operators in terms of vertex operators 
associated with the free boson model. Our construction naturally allows one to treat as well the case of time 
dependent correlations and allows one to recover all the features of an effective field theoretic description which 
are at the very base of the non-linear Luttinger liquid model ifT^ . In this respect, our method, allows one to justify, 
starting from the first principles, the use of the non-linear Luttinger liquid model. 
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A Special functions 


In this section we recall the definition and main properties of the special functions that are of use to our analysis. 
The Euler T-function is defined as 


r(z) = J'^ ■ dt and safisfies r(z -i- 1) = zr(z) . 
0 


(A.l) 


The Barnes G-funcfion is a generalisation of fhe T funclion, in fhe sense fhaf if safisfies fhe functional relation 
G{z + 1) - r(z)G(z). If admits the integral representation 


G{l+z) = (V^-r(z))" • exp 


z(l - z) 


Z 

-f 


lnr(5)-d5 


(A.2) 


from which one can deduce the reflection formula 


G(l-z) 

G(l+z) 


z 

/ 


= (In) ^ ■ exp < I ns cot(7r5) • di > . 


(A.3) 


The latter readily implies the relation 
1 + z, 1 - z - 1 


G 


1 - z, 1 + z + ' 


/sin[7rz]x^ 

[ -) •(-!) 2 

\ 71 / 


(A.4) 


Above we have introduced the so-called hypergeometric tike notation which will be used for products of ratios 
of Euler T or Barnes G-function, e.g. 


{v.}" 

{Wa}^ 


= T 




n r(Vfl) 

a=\ 

m 

n r(H'fl) 
<3=1 


and G 


{Van 


n G(Va) 

<3=1 

m 

n G{Wa) 

<3=1 


(A.5) 


B Some integrals of interest 

In this appendix, we compute some integrals that are of direct application to our study. Namely, given p,h,t € N*, 
define 




r dZdr (1 - gj/Z) 

J (linn Z-T (1 -t/co) 

ld>N>M 


(B.l) 


and 


(B.2) 


r(2)r I 'i - X t ^ (1 - r/m)^ 

M>ld>M 

Two ofher infegrals appear in fhe course of fhe analysis, namely 

7(1)/ I X _ X T~^^P-i (1 - T/gjy ~(2). I , _ ^ . V- 

J(2i7r)2 {-T (l-coizn 9 (linn (1 - m/rf ’ 


ld>l"l>lh 


^d^dr r/’-iX-Ul - 
M>id>i9 
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These two new integrals are however related to the first two introduced above. Indeed, one has 

and jg(y | w) = \ oj) . (B.4) 

The first equality is straightforward whereas the second follows from the change of variables (r, i-> {a/lx, o/ly). 

Lemma B.l The double integrals I w) and I^^piy \ oj) can be explicitly computed as 


hh-1") = 


h-t 


sinfTTv] • OJ' 
n{t — h + v) 


h-V, t + V 

h , t 


(B.5) 


and 


J!J(v I m) 


sinfTTv] ■ 0 )^ ^ P t h + v, p-v \ 
n{h + p - \) \ h , p ) 


(B.6) 


Proof — 

The calculation can be done by means of the method proposed in Appendix A of |T| which allows for an 
effective separation of the integrals. First of all, one observes that 


+ T5r){-^ • C 


(l-rr 




Kv-1 


= vr 


(i-ro 

(1 - 


(B.7) 


Then, upon a rescaling of the integration variables, one gets 




h-t 


dr r ^ 


• ^(1 /iT } X . 


h — t — V J 2\n^ 2i7rT 1 ^ - r (1 - r) 

id>i>M 


(B.8) 


An integration by parts then leads to 

,h-t 


I Oi) = 


-OJ 




rd^dr ,_ 2 (i-ro^ 

h — t — v J (2m)^ (1 - 

ld>IM>|r| 

h-t 


-OJ 


h — t — V 


V-I d^ 


2i7r 


dr 


rxi - X (t) T-\\ - T)-^-^ ■ — . (B.9) 


21^- 


id<i 


|r|<l 


The two integrals can then be evaluated by means of the series expansion 
sinfTTcr] 

«>0 


(1-zf = 


n 


7 , 


a +\ , n — a 
n + \ 


(B.IO) 


hence leading to (IB.51) . 

The integral I^^hv \ oj) is computed along the same lines on the basis of the identity 


+ 'Pdr)[j^ 


y^-T {\-cy 

The details are left to the reader. 


(1-t)7 (1-t) 

- VT 


v-1 


(1 - 


,y+l 


(B.ll) 
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C Proof of Proposition 11.11 

The exponents of current operators preserve the charge of a state. Therefore, the expectation value 

{Jnp-n, ) is non-zero only if 


Up - rih =nk-nt . 


(C.l) 


This explains the occurrence of the Kronecker symbol in (11.361) . 

Observe that the role played by and in the whole construction of the Hilbert space h is symmetric. In 
particular, this ensures that the result will be invariant under the transformation i/^p i-> We 

stress that this transformation leaves the vacuum 10 ) unchanged. Since this transform maps the local current 
operators as 7^. -7^, we get the symmetry 


( T <T \ - r-1 T |p-./-(’'’")p-/'+(''’")| <T \ 

\*Jnp\nh\^ ^ \*Jn}i\nt/ — \ \o/«/,;noF ^ \*Jn!\nk/ 


{C.l) 


where the sign prefactor issues from the change in the order of products of the t^’s and i/r*’s. It is easy to check 
that indeed, as soon as the condition rip - ni, = - rit is imposed, the rhs of (11.351) is indeed invariant under the 

simultaneous transformations 


{iTnp'.ni,^ ^ {^ni,',npi iTnt.rif ' 


(C.3) 


followed by a multiplication by (-1)”'+”''. The above thus ensures that it is enough to establish the validity of the 
representation (11.351) solely in the case when rip - rih = rik - rit > 0. Hence, from now on, we focus on this case 
only and, so as to lighten the notation, we introduce the shorthand notations 


tifi = n rip = n + { rit = d and = d + £ . (C.4) 

In order to compute the expectation value, we first reduce the problem to a simpler one where the application of 
Wick’s theorem is straightforward. This boils down the problem to the computation of determinants whose entries 
are given by hi- and quadri-linear expectation values in the fermions. The expectation values can be expressed in 
therms of the integrals computed in Appendix I bI Building on these results, we compute the various determinants 
by using their relation to Cauchy determinants. 


• Reduction to a simpler problem 

Observe that, for any r > 0, one has 

= <01 . (C.5) 

Let ri,..., be distinct integers satisfying to the constraints 

ry > max{/ia : a = l,...n} for any b = l,...,£. (C.6) 

Then, in virtue of (IC.5I) . we get 


{Jnp-,n„ ) - <0| f] ) 


<3=1 


= (-1)^”( 0 • if'-K • • • <A-/n f-n ''' ) . (c.7) 


Pn-¥i 
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We now move the operators il/% through the exponents of current operators. For this purpose, we first represent 
in terms of the conjugate field by means of a contour integral around 0 and then we apply the exchange 
relation (11.251) : 


(Air 


(1 - T/m) 




= . (J)* ... (J,* O* 


<1-Wrr 


Thus, upon extending the sequence ha as h„+i, - rf,ioxb - 1,..., A and after setting 

K = 'P-ta a = l,...,d ^ a=l,...,£ and ^ iAa„-i a = l,...,d + e, (C.9) 

the expectation value we started with reduces to 

< Jn,-n, |e^-(''’-)e) - (-1)^(”+'^)( w*, ■ ■ ■ m 10) 

= 0 ) . dew [Mab] . (C.IO) 

The set J'n+e;n+{ appearing above is defined in terms of the extended sequence {haT^^ as 

JnM = [{Pa]T' ; iha}T'} ■ (C-H) 

We also specify that we have obtained the second line in (1C. 101) by applying Wick’s theorem. Finally, the matrix 
M arising in (1C.101) reads 

< |e^-(^’-)e J10 ) 


^ab ~ 




for a = I,... ,d and 


Ma+db = (Dt' 


,„_i(l-T/mr < Jn+f;n+b 10 ) dr 


(1 - oj/ry iJn+f-n+t |e A(v.^)| Q ) 


2i7r 


(C.12) 


(C.13) 


for a = 1,..., A. 


• Evaluation of {J'n+t.n+t 0) 

The expectation value {J'n+(;n+e 0 ) can be computed by applying, again, Wick’s theorem: 

< 0) - det„+^ [< 0 0)] (C. 14) 

where we have used that e‘^+‘^’'’"^|0) = |0) and (0|e-^= (0|. The entries of the matrix arising in the 
determinant coincide with one of the integrals studied in Appendix |B] Indeed, one has 


<0|.A;-i<A-a • 0 } = ^ |T'*(t)T'(^) • e^-(’'’->|0) 


|r|>iq 


drd^ 

(2i7r)2 


|r|>ld>|t^| 


(1 - gj/^y dTd{ 

T-{ (1 - oj/tY ( 2 i 7 T) 


,2 = I • <C. 15 ) 
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In the intermediate calculations, we have used the exchange relation (11.251) along with 


<oi'F*(Tm^)io) - 




for |t| > 1^1 . 


Hence, we arrive to the determinant representation 

- dew I m)] . 

• Evaluation of {Jn+e-n+e 0) 


(C.16) 


(C.17) 


Since the operator e» 




is a group-like element, Wick’s theorem ensures that 


< <A^-i 10) 


- (0|e^-('’’-Vl*<A^-i|0)-deW^ 




(0|eA(v.^Vl,Wi|0) 

The two-fermion expectation value can be computed analogously to (IC.15I) 

drd^ 


mT\>m 

The four point function is recast as 


^-T (l-^/mf (2i;r)2 


= l?\y I ■ 


(C.18) 


(C.19) 


( 0 WU Wi 10 ) 




jt+\^k ("I _ ailzYi^i - Tliof 
|z|>|f|>|w|>|r|>|f| 




dzd^'drd^ 

(2i/r)4 


. (C.20) 


A straightforward computation shows that 

<oiT'*am^)'p*(rm^)io) - 


T • Z 


T • Z 


(^ - t) • (z - (r - f) • (z - 0 


(C.21) 


provided that the variables satisfy to the ordering |z| > l^"! > |t| > |^|. After inserting the above expectation value 
in (lOOl) . we are led to 

< 0 |^;_i , Wi10 ) = Ji;\v I m) • (v | m) + lf{v \ ca) ■ lf^{v \ co) . (C.22) 

All this allows us to represent the expectation value of interest as 


< J'«+r;n+r >l'k-i\0) = lf{v \ m)• det„+^ (v I m) + - 


if{y I a;) 


(C.23) 


"'(2) 

This representation involves the determinant of the matrix Nab = ^ p I <^) perturbed by a rank 1 matrix. Since 
the matrix N is of Cauchy type, its inverse is explicitly computable and takes the form 


(W-) 


nu) 


i-Pb-K 


Pb 


hf. 


n+€ 

1 p=i 


n+t 

n {(K + pfo - 1) • {ha + Pc - 1)} 


ah sin[7rv] \pfo + V, ha-V ) ha+Pb-l Ucl^iPc - Pb) ' UcY{(>^c - ha) 

i^b +a 


. (C.24) 
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Thus, recalling the representation (IC.17I) we are led to 

< «A^-11 0 ) - lf{v I CU) • ( |e A(v,-)e 0 ) • Sr,k{jn^t,n^i) . (C.25) 

The set function St,k{Sfn+t-,n+?j takes into account the contribution issuing from the rank 1 perturbation of N and 
takes the form: 


'St,k{<T>i+i\n+t^ — 1 + ^ 




a,b=\ 


-rliVlw) 


n+i 

= >-z 


t + k - l 


n+l , . 

n [{he + Pa - 1) • {hb + Pc - 1)} 

C=1 ^ 


a A -^>’ + V){k - Pa -V) YlTJliPc - Pa) ' Y\T=\{hc - hb) 

+a i=b 


= 1. 


t "¥ h — 1 


J T\n j T\n (^ + z - 1)(? - Z + y){k - ^ - v) 
r(ip„i) r(|A„i) 


n+t ^ 

ni 

c=l 


{he + s - l) ■ {z + Pc - 1) 

{Pc - s) ■ {he - z) 


. (C.26) 


Above, T{{pa}), resp. r({/ja)), stands for a small counterclockwise loop around p\,... ,Pn+e, resp. hi,... ,h„+f, 
that avoids all other singularities of the integrand. Note that the contour integral should be understood as an 
encased integral. As a consequence, the z-integrand does not contain poles at z = 1-5. The integral can be 
evaluated by taking the residues located outside of the contours of integration: 


St,k{Sfn+i-,n+£^ 


= 1 + 


d5 


t k — 1 


lin (^ + ? + y - l){k - s - y) 


r(lpd) 


n( 

n+^ 


ni 


{he + S - \) ■ {1+ V + Pc - l) 

{Pc - S)- {he-t-v) 

{he + k — V — \) ■ {t + V + Pe — V) 


{Pe - k + V) ■ {he - t - V) 


(C.27) 


All-in-all, we are thus led to the representation 

< |e A(v.^)e 1 0 ) _ ^( 2 ) 




h?(v I <.>). n 


C=1 


(/jc -I- k - y - 1) • (t -I- V -I- Pc ~ 1) 

{Pe - k + V) ■ {he - t - V) 


(C.28) 


• Evaluation of Ma+db for a = 

It follows from Wick’s theorem that one has the representation 

{JnU-nU |e^-(''’'^)eA(v.-)|0) . Ma^db 


(2)/ % 

(i-T/mr4w 


-• det„+/ 

(1 - m/rr 




dr 

(C.29) 

2i7r 


where q, m label the entries of the matrix whose determinant is to be evaluated while the function gf\T) is defined 
by means of a contour integral: 


Sk 


( 2 ) 


(T) = <0|e=^-(''’'’^V(T)<A^,_i|0) = 


N>M>ld 


T-r^ {\-acoy 

T - ^ (1 - tIojY lin 


(C.30) 
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(C.31) 


The expectation value in the determinant is computed along the lines already discussed at length and reads: 

< 0 |^;_i 10 ) - I m) • + J^Cv | m) • gf (r) 

where the function is defined as 


g^\T) ^ 






{I-UJIO^ d^ 
(1 - tIojY lin 


(C.32) 


ld>M>|r| 

Upon factorising the matrix {N)ab - I from the determinant in (IC.29I ) we are led to the representation 


Ma+db 


■c 


-h 


|r|>IM>ld 


T-i 


(1 - qcjy d^dr 

(1 - ojItY ' (2i7r)2 


+ 


n-Yl 


m,q=\ 


. j(l) 

Pqh 


(y I m) 





{l-oj/0^ d^dr 
(1 - oj/tY (2mY 


(C.33) 


The first integral vanishes as can be seen by deforming the r-contour to oo. Likewise, by deforming the t contour 
to oo in the integral arising under the double summation symbol, one solely gets the contribution of the pole at 
T = ^. The latter yields the Kronecker symbol 6m,n+a which is there so as to enforce the constraint that m should 
be such that h,„ = ra- After some algebra, one recasts Ma+db in the form 


Ma+db = ■ r( 

V hn+a - y 

where 


kb -V 
kb 


1 ) 


n+€ 

0 Q^n-k-a "I" Pc 

L-1_ 

ntil - hn+a) 

i^n+a 




n-¥a,b 


(C.34) 


n-Yl 

n+{ . 0 (^c Pq ~ 1) 

^ ^ V_^_^2^_ 

^ fe - Pq-y)- (hn+a + Pq - 1) n”=[(Pc “ Pq) 

*q 


r -1 PI (/ic + - 1) 

di 

, n+£ 

n 

J {kb- S-v)- {hn+a + 5-1) {Pc - s) 

2i7r 

hn+a +kb-y-l 1^1 


r(lpd) 


{hc + kb-y- 1) 
{pc -kb + y) 


(C.35) 


• Calculation of deU+f [Ma b\ 

Summarising the results obtained so far, we get 

_ sinfTTv] • ! kb-y, U + V 

~ n{ta+kb-\) \ kb , ta 


for a = I,.. .,d and b = I,... ,d + £ while 


n+e 

n( 


{Pc + tg + y - 1) ■ {he + kb - y - 1) 
{he -ta- y){pc -kb + y) 


Mg+d b — 




V 1 \ ^n+a ^ •> 


hyt 


kb-y \ n”=f [{he + kb-y- l){hn+a + Pc - 1)1 
' ' Yret{{Pe -kb + y)- ntil ^he - hg+a) 


(C.36) 


(C.37) 
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We now recast the matrix Mat into a uniform representation. Having this in mind, we introduce two sequences 
a a, and with a = \,... ,d + €. The sequences a a reads 


(Xa 


sin[7rv] . p/ + V \ |-4 r (Pc + tg + v- 1) ] 
7T ■ 0)^0 I ta I 11 I {he - ta-v) / 


and 

’-n+s Pc ~ 1) 

{he — /jg+i) 


ad+s ^ -OJ 


•ri 


K+s - V ) YVel] 


where the index a runs over {1,..., <i) while the index 5 runs over . .,i}. Furthermore, let 


Pb - 


i • fl I + h - y - ^) \ 

V h / I {Pe -kb + v) ) 


for b = I,... ,d + £ . 


(C.38) 


(C.39) 


(C.40) 


Then, upon extending the sequence for a = r/ + 1,..., J + ^, as 


ta+d 


~ ^n+Q y 


(C.41) 


we see that M is closely related to a Cauchy matrix 

Mat = , , for a,b e . ,d + £} . 

~ ^ 


(C.42) 


• Synthesis of the calculations 

By gathering all together the results obtained in the course of the proof, one arrives at the representation 




b_ubn+d) I - sin[7rv] • ^ ( h - v, Pa + v \ 

^ ^ [ n{pa + hb-l) \ hb , Pa iJ 

d+£ 

X f] 1». ■ Al ■ [pT^l ■ [^^TfeTrl ■ (C-«) 


It is then a matter of straightforward calculations to obtain the representation (I1.35I) . In particular, all the depen¬ 
dence on the auxiliary integers {ra}^^ completely disappears. ■ 


D The multi-point restricted sums 

The multi-point restricted sums were first introduced in 1(2711 . This name refers to a multi-dimensional summation 
formula that is necessary for the computation of the localised form factor series expansion. More precisely, given 
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kil > • • • > \Zr\ and generic complex numbers vi,..., it holds 


r-1 +00 r—1 r—1 

y{{VaY{,{Za}\) = Y[[ Z Z 

s = \ „(■!) „(s)_A ,T (A .5=1 ' ^ S=2 


ny,nf=0 J W^„W 


r-1 


i=l 


^ nK^) 


2 


G 


1 + ~ Vi, 1 + /’i + Vi+l 


1 - Vi, 1 + Vi+1 

xfi^C! 


s=2 


l+Vijl+Z^i—^ — ^i"t Vi 
+ Vi, 1 + £s-i ~ Vi 


no-S)' 


Zb \(ya+Ka){Vb+Kb) 


■ (D.l) 


b>a 


The function K appearing above takes the form 


K(jn,-n„ I V,?7;z) 


sin[7rv] 

n 


sin[7r77] \”* 
n ) 


rip rib 

n (Pa - Pbf ■ Yliha- hf 

a<b a<b 


rip rib 

n t\{Pa + hb- 1)2 
a=\ b=\ 



Pa-V,Pa+P 
Pa, Pa 



K + V,ha' 
f^a, f^a 


This summation formula has been argued in GTI I on the basis of comparing the l&rge-N asymptotic expansion 
of a Toeplitz determinant generated by a symbol having Fisher-Hartwig singularities obtained by two different 
means. In this section of the appendix, we derive the above formula by using the free fermion formalism that we 
have developed. For this purpose we compute the expectation value 


{0\y{vi,Ki I Z\)---'y{Vr,Kr \ Zr)|0) 

in two ways. First of all, by inserting the decomposition of the identity 


(D.2) 


2] 2] \YTnp-,nb){Jnp-,nb \ = id (D.3) 

rip,rib 

between each vertex operator and this, for the specific choice of the sequence Ka ■ 

Ks - (s-i - £s where £q - £p - 0 . (D-4) 

Second, by applying repeatedly the exchange relation 

g_/+(v,w). - |l - —J'’ ■ . (D.5) 

satisfied by fhe current operators. The second method yields 

{0\y{vuKx\Zl)---y{Vr,Kr\Zr)\0) = (D_6) 

b>a 
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while the first one leads to 


r 

{0\y{n,K,\Zx)---y{Vr,Kr\Zr)\0) = Z H {< 


„«-„«= JW. (.) ^=1 


, (5+1) (5+1) 

Ks+i+n'p -n\ 


Above, we have adopted the convention that 


(0) (0) (r) (r) 

np> = n^> = ny = nl> 


= 0 viz. \J„( 0 ).„( 0 )) = \ J^(r).^(r)) ^ |0). 

ftp i'll, "p 


(D.8) 


The constraints on the differences in fhe numbers of parficles and holes in fhe excifafions of fhe infermediary sfafes 
issue from charge conservafion requiremenfs and can be recasf, due fo fhe paramefrisafion (ID.4I) . 


n^p = (s 5 = l,...,r - 1 . 

Then a slraighlforward calculafion recasfs fhe sum as 

{0\y{vi,Ki\z\)---''V{Vr,Kr\Zr)\^)-C{[VaY{,{ZaYi) = yiiVaYv iZaYi) 
where fhe consfanf C{{vaY{, {Zalp takes the form 


(D.9) 


(D.IO) 




.V=I 


\ks(Ks-1)+Ks{s 


^5-1 ^^5-1'*'^^ I ^5(*^5 + l) 


(- 1 ) 


(D.ll) 


By using the reduction property of the Barnes function (IA.4I) and the parametrisation of Ks in terms of £s, one gets 
that 


5=1 ' ' 


sin[7rvj ^-4-1 _ 


^5— 1 (^5- 1 ^ 1 Xs (^5 + 1)1 

2 2 


r-1 

n 

.v=l 


n 


G 


\ + is — Vs,^ + is Vi+l 
1 - y^, 1 + y,+i 


r-l 


s=2 


1 + y^, 1 + is-\ ~ + y^ 

+ y^, 1 + is-\ ~ Vs 


, (D.12) 


and 


r r-l 

yKs(Ks-l)+Kst: 


ni 


FiHx)'") = n{(^) 


.V=l .V=l " i=l 

Putting all of this together leads to (ID.ll) . 


tslhrl) 

2 


(D.13) 
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